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ABSTRACT 

In this paper we study the oxygen chemistry in the C-rich circumstellar shells of IRC+10216. The 
recent discoveries of oxygen bearing species (water, hydroxyl radical and formaldehyde) toward this 
source challenge our current understanding of the chemistry in C-rich circumstellar envelopes. The 
presence of icy comets surrounding the star or catalysis on iron grain surfaces have been invoked to 
explain the presence of such unexpected species. This detailed study aims at evaluating the chances of 
producing O-bearing species in the C-rich circumstellar envelope only by gas phase chemical reactions. 
For the inner hot envelope, it is shown that although most of the oxygen is locked in CO near the 
photosphere (as expected for a C/O ratio greater than 1), some stellar radii far away species such as 
H2O and CO2 have large abundances under the assumption of thermochemical equilibrium. It is also 
shown how non-LTE chemistry makes very difficult the CO— >H20, CO2 transformation predicted in 
LTE. Concerning the chemistry in the outer and colder envelope, we show that formaldehyde can be 
formed through gas phase reactions. However, in order to form water vapor it is necessary to include a 
radiative association between atomic oxygen and molecular hydrogen with a quite high rate constant. 
The chemical models explain the presence of HCO + and predict the existence of SO and H2CS (which 
has been detected in a A 3 mm line survey to be published). We have modeled the line profiles of 
H2CO, H2O, HCO + , SO and H2CS using a non-local radiative transfer model and the abundance 
profiles predicted by our chemical model. The results have been compared to the observations and 
discussed. 

Subject headings: astrochemistry — circumstellar matter — molecular processes — stars: AGB and 
post-AGB — stars: individual(IRC+10216) 



1. INTRODUCTION 

IRC+10216 is a low mass AGB star losing mass at a 
rate of 2-4xl0 -5 M Q yr _1 in the form of a molecular 
and dusty wind that produces an extended circumstellar 
envelope (CSE). The processes of dredge-up that occurs 
during this evolutionary late stage alter the elemental 
composition in the stellar surface (that is roughly solar, 
C/0<1, during the main sequence phase) resulting, in 
the case of IRC+10216, in a C/O ratio greater than 1, 
in what is known as a carbon star. 

The physical conditions in the vicinity of the photo- 
sphere of an AGB star (typical temperatures of ~2500 
K and densities of ~10 14 cm~ 3 ) make the material to 
be mainly molecular, with a composition determined by 
local thermodynamic equilibrium (LTE). It i s now well 
established since the pioneering work of frsiil lfl973T) that 
the C/O ratio completely determines the kind of chem- 
istry taking place. The high stability of CO makes this 
molecule to have a large abundance locking almost all 
of the limiting reactant and allowing for the reactant 
in excess to form either carbon bearing molecules when 
C/0>1, and the opposite, i.e., oxygen bearing molecules, 
when C/0<1. This has been extensively confirmed by 
observations. A look to the list of molecules detected in 
IRC+10216 confirms that t hey are mostly carbon b ear- 
ing species (see table 1 in iCernicharo et al.l 120001 and 
§3). The only oxygen bearing molecule in C-rich AGB 
stars with a signific ant abundance, apart from CO, is SiO 
ijMorris et alJll975j) . 

However, some other O-bearing molecules such as 
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HoO HMelnjck eT^H I200T1 EWawa et al.l EflOrl. OH 
l|Ford et al.H2003|) . and H^CO l|Ford et alJl2004|) have re- 
cently been detected in IRC+10216. The existence of 
such molecules has been interpreted as the result of evap- 
oration of c ometary ices from a K uiper belt analog in 
IRC+10216 i|Ford fc Neufeldl 12001]) . The luminosity in- 
crease of the star, intrinsic to the red giant phase, would 
have caused the ice sublimation and subsequent release of 
water vapor to the gas phase in the circumstellar shells. 
OH would be produced when water is photodissociated 
by the interstellar UV field in the unshielded outer enve- 
lope and H2CO would be the photodissociation product 
of an unknown parent molecule produced by sublimation 
of the ice mantles of these comets. 

An a l ternat ive explanation has been proposed by 
iWillacvl l)2004|) in which H2O would be produced through 
Fischer- Tropsch catalysis on the surface of iron grains, 
that would be present in the expanding envelope due to 
condensation from the gas phase of some fraction of the 
available iron. Fischer- Tropsch catalysis breaks the CO 
bond and produces H 2 and hydrocarbons such as CH4. 

O-bearing species have also been detected in the C-rich 
proto planetary nebula CRL618 l)Herpin fc Cernicharo! 
2000). However, th e production of th ese species has 
been interpreted by ICernicharo! (|2QQ4l) as the result 
of a rich photochemistry in a region of high den- 
sity (~10 7 cm -3 ) and temperature (200-300 K) where 
large complex carbon-r ich molecules are also produced 
l|Cernicharo et alJl200ll) . 

The aim of this paper is to investigate whether oxygen 
bearing molecules could be produced in the carbon-rich 
expanding gas of AGB stars by non-LTE mechanisms. 
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LTE calculations provide a good estimation of molecular 
abundances in the vicinity of the photosphere but as the 
gas expands the temperature and density decrease sig- 
nificantly and the chemical timescale increases, making 
chemical kinetics dominant in determining the molecular 
abundances. We describe the model of the circumstellar 
envelope in §2. The reaction network is discussed in §3. 
The results of the chemical model are presented in §4 
together with comparisons of radiative transfer calcula- 
tions with available observations. The conclusions are 
given in §5. 

2. THE CIRCUMSTELLAR ENVELOPE 

We assume spherical symmetry for the circumstellar 
envelope. In order to calculate the molecular abundances 
at different radii in the CSE, we follow the history of a 
volume element of gas with a given chemical composi- 
tion travelling from the photosphere (r=i?*) to the end 
of the envelope (r^lO 18 cm). We build a system of dif- 
ferential equations which integration yields the temporal 
evolution of the density of each gas species. The time 
has to be interpreted as radial position provided the gas 
travels outwards at a given velocity. Different processes 
are considered during this travel depending on the posi- 
tion in the CSE. Kinetic temperature (which determines 
the reaction rate constants) and total gas density radial 
profiles are needed to solve the system of equations. The 
values adopted for the different parameters used to model 
the CSE are given in tabled 

For model purposes the CSE is considered to consist of 
three different regions: the innermost region, the inner 
envelope, and the intermediate and outer envelope. 

2.1. The innermost region 

This zone corresponds to the region between the pho- 
tosphere and ro (see Fig. QJi) . The adopted physical 
conditions are: 

- The temperature is considered to vary as a power law 
of the radius r 



T(r) = T* x I — 



(1) 



where T* is the temperature at the photosphere and i?* 
is the star radius. 

- The density profile is given by hydrostatic equilib- 
rium. Considering the above temperature law, the gas 
density can be expressed as follows in terms of the value 
at r , n(r ) 



n(r) = n(ro 



exp 



GM*rriH(J< 
'kT,R?{l-a) 



.(«-!) „(«-!) 
' 

where k is the Boltzmann constant, G is the gravitational 
constant, M* the mass of the star, to# is the mass of a 
hydrogen atom and fi the mean molecular weight of gas. 

The high densities and temperatures in this region al- 
low to assume that molecular abundances are given by a 
LTE calculation. 

2.2. The inner envelope 

The gas travels from ro to r c (see Fig. QJi), where 
r c stands for condensation radius. At this distance, 
grain formation is supposed to occur. Radiation pres- 
sure on grains, together with momentum coupling of 
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Fig. 1. — (a) Three different density profiles considered for 
the inner and innermost envelope: — MH [K n =l, n(ro)=3.7x 10 14 
cm" 3 ], -MI [Jf n =1.8565, n(r )=3.7x!0 14 cm" 3 ] and -ML 



[K n =l, ra(r )=1.2 X 10 1 



all them use a r law for r>r c 



(b), (c) and (d) show the abundances of H and H2, relative to 
total number of hydrogen nuclei [n(H)+2n(H2)], as given by LTE 
(dotted lines) and by chemical kinetics (solid lines) for the three 
density profiles: (b)— ►MH, (c)— >MI and (c)— >ML. Note the strong 
dependence of the H/H2 ratio on the densities considered. 

gas and grains, make the gas accelerate up to a termi- 
nal expansion velocity which remains constant beyond 
r c . The mechanism responsible for the transport of the 
gas from ro to r r is somewhat controversial. Hydrody- 
namical models l)Bowenl|l988l) have shown that pulsa- 
tion s,! driven shocks can gradua lly move the gas up to 
r c . iWillacv fc Cherchneffl l)1998fl have applied this ap- 
proach for modelling the chemistry in the inner envelope. 
An alternative mechanism could be that the gas close to 
the photosphere pulsates around an equilibrium position 
with an associated radial velocity very similar to the es- 
cape velocity. This could lead to a scenario in which 
shells of gas that are in levitation can eventually escape 
from the surroundings of the photosphere and some of 
them could reach r c . 

- The temperature profile in this region is given by 
equation 1. 

- For the density profile we follow the treatment 
of a shocked extende d region of a C-rich CSE of 
ICherchneff eiTaH (Jl992) 



n(r) = n(ro)exp< —K, : 



GM»m g /i(l - 7 2 
kT*R«{l-a) 



(3) 

where 7 is defined in ICherchneff et al.l l)1992fl and K n is 
explained in the next paragraph. 
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The reference value n(ro) is not evident neither from 
observations nor from models. Previous models of the 
inner envelope of IRC+102 16 ijCherchneff et alJ 119921 
iWillacv fe Cherchnefll [l998) have used high values of 
n(ro) based on hydrodynamical models for periodically 
shocked Mira-like stars, which have density profiles as 
an output. However, with such high n(ro) values, the 
extrapolation to larger radii produces densities for the 
outer envelope much higher than those obtained from 
the law of conservation of mass, in which the absolute 
density is well determined by the mass loss rate and the 
expansion velocity (see equations 5 and 6 in §2.3). To 
reconcile the outer envelope density law with the high 
values in the inner regions we introduce the arbitrary 
factor K n in equation 3. We consider three different den- 
sity laws for the inner envelope, shown in Fig^,, whose 
parameters n(ro) and K n are given in table MI (I 
for intermediate density); MH (H for high density), also 
considere d for comparing our results with the previous 
model of Will acv fc Cherchnefll l)1998t) although it over- 
estimates the density for large radii; and finally ML (L 
for low density). We point out that the chemical time 
scale strongly depends on the density and also the H/H2 
ratio is completely different in a low or high density sit- 
uation (see Figs. H 3 ,^ and^i). We will comment this 
in §4.1. 

Another issue is how the gas travels from r$ to r c . As 
discussed above, different mechanisms could be operat- 
ing. The simplest assumption to model the chemistry is 
to consider that the gas travels at constant velocity in 
this inner envelope. When using this approach we will 
assume a velocity of 1 km/s. 

- SiO condensation onto SiC grains. Considering 
non-LTE effects, we investigated whether deposition of 
Si species (that are specially refractory) on SiC grains 
could affect the oxygen chemistry because SiO is the 
second more abundant O-bearing molecule. The scheme 
for grain formation in IRC+10216 could be as follows. 
The major type of grains are carbonaceous, most proba- 
bly amorphou s carbon (A.C.) as in dicated by fitting the 
IR spectrum ijBaenulo et al.lll995j) . These grains would 
originate in condensation processes involving C2H2 (the 
most abundant C-bearing molecule after CO) that 
would take place at ^1000 K (^4 R»). A minor type 
of grains would be SiC, as indicated by the 11.3 nm 
featu re in the IR spectrum l)Lorenz-Martins fc Lefevrel 
1993), that would condense at ~1500 K (~2 R»). This 
condensation sequence is due to an inverse greenho use 
effect in both types of grai ns (see IMcCabejf T982 ] for 
a theoretical approach and Frenkl ach et alJ 119891 for 
experimental evidence). Although a more correct 
approach would be a chemical mechanism leading to 
the formation of grains from gaseous molecules, for our 
purposes we describe the condensation of SiC grains 
(A.C. grains are not considered because we are just 
interested in deposition of Si species on SiC grains) as a 
two-step process: (1) nucleation from SiC gas molecules 
in the context of homogeneous nucleation theory, which 
is basically described by the nucleation rate J» (number 
of condensation nuclei formed per unit time and unit 
volume); and (2) growth of grains by accretion of gas 
phase molecules, that is described by the grain radius 
a gr . The calculation of these two magnitudes is detailed 
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Fig. 2. — Parameters for a volume element of gas that suffers 
a history of 5 shocks: (a) trajectory; (b) temperature history; (c) 
density histories for models MH, MI and ML. (d), (e) and (f) show 
the H and H2 abundances, relative to total number of hydrogen 
nuclei, calculated by chemical kinetics for the three density laws 
shown in (c). 

in Appendix A. 

- Shocks. An alternative to the mechanism of the gas 
travelling at constant velocity in the inner envelope is 
to consider the effect of gas driven by shocks. Molec- 
ular abundances are considerably affected by the steep 
changes in temperature and den sity associated to the 
shocks . We follow the approach of Willac v"fc Cherchnefll 
( 1998) but in a more simplistic way. The history of a 
volume element of gas to reach r c from ro after a certain 
number of shocks is represented in Figs. |2t-, I2t> and[2t, 
and is described below . With one shock and subsequent 
relaxation, the gas moves from a position r% to r 2 (r 2 >ri) 
following a harmonic trajectory. The gas in r\ with tem- 
perature T\ and density n± suffers a shock and we assume 
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TABLE 1 
Model parameters 



Abundances relative to H2 
for parent species entering 
IRC+10216 parameters Neutral molecules considered into the outer envelope 



[1] 




6.5 10 13 


CH 
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N 2 


SiH 




H 2 
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M,/M Q 


2 


CH 2 


O2 


NH 


SiH 2 


[18] 


H 


4.5 10 


-3 


[2] 


ro/R* 


1.2 


CH 3 


OH 


NH 2 


SiH 3 


[10] 


He 


2.0 10 


-1 


[1] 


r c /R, 


5 


CH 4 


H 2 


NH 3 


SiH 4 


[7] 


CO 


6.0 10 


-4 


[1] 


T»/K 


2320 


c 2 


H0 2 


CN 


Si 2 


[19] 


C 2 H 2 


8.0 10 


-5 


[3] 




0.6 


C 2 H 


H 2 2 


HON 


SiC 


[19] 


HCN 


4.9 10 


-5 


[2] 


n(ro) MH ' MI /cm- 3 


3.7 10 14 


C 2 H 2 


CO 


NS 


SiCH 


[10] 


N 2 


8.0 10 


-5 




n(r ) ML /cm- 3 


1.2 10 12 


C 2 H 3 


co 2 


SH 


SiCH 2 


[13] 


C 2 H 4 


2.0 10 


-8 




t/MH.ML 


1 


C 2 H 4 


HCO 


H 2 S 


SiC 2 


[11] 


CH 4 


3.5 10 


-6 




K MI 

n 


1.8565 


C 2 H 5 


H 2 CO 


cs 


Si 2 C 


[11] 


NH 3 


1.7 10 


-7 


W 


P 1 days 


650 


C 2 H6 


CH 3 


HCS 


SiN 


[11] 


SiH 4 


2.2 10 


-7 
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7 


0.89 


c 3 


CH 2 OH 


H 2 CS 


SiNH 


[14] 


H 2 S 


6.0 10~ 


10 


[3] 


/i(r<r c )/amu 


1.7 


C 3 H 


CH 3 OH 


S 2 


SiO 


[15] 


CS 


1.2 10 


-7 


[5] 


v exp /(km/s) 


14.5 


C 3 H 2 


HCCO 


CS 2 


SiOa 


[16] 


SiS 


4.3 10 


-6 


[6] 


M/(M yr" 1 ) 


3 10~ 5 


C 4 


NO 


SO 


SiS 


[17] 


SiC 2 


5.0 10 


-8 


[7] 


distance/pc 


150 


C 4 H 


NCO 


so 2 




[11] 


SiO 


8.0 10 


-7 


[8] 


riooo(r- = 10 16 ) 


12.7 


C4H 2 




ocs 
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that the immediately post-shocked gas increases its tem- 
perature by a factor 2 a nd its density by a factor 10 (typ- 
ical values obtained bv lWillacv fc Cherchnefflll998l from 
hydrodynamical considerations). Then the gas relaxes 
its temperature and density exponentially until reaching 
the new radius r-i (with temperature T2 and density n^) 
in a time interval equal to the period of pulsation of the 
star P. In this way the gas progressively moves outwards 
until reaching r c , beyond which an expansion at constant 
velocity v exp takes place. 

In Fig. [5] we plot the trajectory, temperature and the 
three density laws introduced previously. Note that again 
the H/H2 ratio strongly depends on the total gas density 
(see Figs. 01, andEf). 

2.3. The intermediate and outer envelope 

It extends from r c . The gas expands adiabatically at 
constant velocity v exp . The molecules in the very outer 
envelope are no longer shielded by the dust of the CSE 
against the interstellar UV field, and are then photodisso- 
ciated. Cosmic rays also play a role ionizing some species. 

- The temperatur e profile is taken f rom the fit of 
M amon et alJ 1^)88) to the results of IK wan fc Linkd 
(|1982() for the outer envelope of IRC+10216 



T(r) = max 



10,14.6 



9 x 10 16 



(4) 



where r is expressed in cm and p takes values of 0.72 for 
r<9xl0 16 and 0.54 for r>9xl0 16 . 

- The density profile is given by the law of conservation 
of mass: 

M = n{r)niH pvexp^r 2 (5) 
which, with the values given in tableland assuming that 



hydrogen is mostly molecular (/i^2 amu), results in 
3.1 x 10 37 



n(r) 



,-3 



r(cm)' 



(0 



- The photodissociation rate Ti for a molecule i at a 
radius r depends on the UV field at that radius: 



4TrJx(r)a l (X)d\ 



(7) 



91.2 



- 2 S" 1 



where 47rJ\(r) is the UV field at r in photons cm 
nm _1 , o~i{\) is the photodissociation cross section and 
the integral extends from 91.2 nm (Lyman cutoff) to a 
threshold value A; that depends on each molecule. The 
calculation of each Ti{r) involves the knowledge of the 
cross section and the solution of the UV radiative trans- 
fer for calculating J\ at each radius r, which will depend 
in a first approximation on the amount of dust that sur- 
rounds the point r. In practice, photodis sociation rates 
are u sually expressed in the literature ijLe Teuff et alJ 
2000) for plane-parallel geometry as a function of the 
visual extinction A v (in magnitudes) measured along the 
direction normal to the infinite plane 



Ti(A v ) = A ie xp(-CiA v ) 



(8) 



with specific parameters Ai and Ci for each molecule. 
For our purposes we have adopted this simple approach 
just correcting for the geometrical difference from plane- 
parallel to spherical. The procedure used for such correc- 
tion as well as the resulting A v radial profile are detailed 
in appendix B. 

3. THE CHEMISTRY 

We cons ider solar element a l ab undances (oxy- 
gen from lAllende Prieto et alJ 120011 and the rest 
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from iCoxl 12000(1 and assume a C/O ratio of 2. 
The chemistry in IRC+10216 is dominated by car- 
bon bearing molecules. F or ca rb on chain radi- 
cals C „H see iTucker et all J1974D: IThaddeus et al.l 
1119851): iGuelin et all (I1978TI: iCernicharo et a.U (|198fif> : 
S fit, aUII 987411 997MCfimicWo fc CuelirJ (1996). 
For carbon chain radicals C 4M see iWilson et al. (|1971|) : 

IGuelin fc Thaddeusl (1T977T) : IGuelin et all (1199? ) For 

cyanopolyynes HC271.+1N s eel Morris et all (119711 119 76): 
iWinnewisser fc Walmslevl l|1978|) : iBell et a I il!982l)~r 
sulfur-carbon chain mo l ecules C„,S see Wilson et alJ 
lfl97l|k ICernicharo et all (fl987^ : iBell et alJ (I1993D. And 
for silic on bearing molec ules SiC „ see TCernfc haro et ahl 
(IT989I): IThaddeus et alJ (fl98^) : lApponi et all (|1999|) ; 
Ohi shi et all l|1989|) . However, since we are just inter- 
ested in oxygen chemistry, in our model we don't include 
very long carbon chains but only the 65 neutral molecules 
given m tabled 

As initial abundances for the inner envelope we take 
those given by LTE at ro, while for the outer shells we 
adopt initial abundances for some parent molecules taken 
from the literature or from a LTE calculation when no 
data are available (see table P). This procedure is mo- 
tivated by the processes occurring in the intermediate 
envelope (specially on grain surfaces) that presumably 
alter the molecular abundances (hydrides such as CH4, 
NH 3 , SiH4 and H 2 S are supposed to form via such pro- 
cesses) . 

Different types of reactions dominate the chemistry de- 
pending on the region of the CSE considered. The tem- 
peratures in the inner envelope are not high enough to 
make ions abundant but allow for reactions with activa- 
tions energies of up to a few tens of thousands of K to 
occur. Therefore, important processes are termolecular 
reactions, its reverse (thermal dissociation) and bimolec- 
ular reactions between neutrals. The outer envelope 
chemistry is dominated by photodissociations, radical- 
molecule reactions without activation energy and radia- 
tive associations. Although all the species detected in 
IRC+10216, except HCO + , are neutral, ionic chemistry 
is necessary to explain the formation of some species. 
Therefore, when modelling the chemistry in the outer en- 
velope we add 80 ionic species (mainly the positive ions 
and protonated species of the neutral molecules of ta- 
blc^) and include the subsequent reactions in which they 
are involved. The rate constants have been taken from 
different sources: databases such as NIS T Chemical Ki- 
netic s Database 1 , UMIST Data base 1999 flLe Teuff et alJ 
2000) and osu.2003 Database l|Smith et alJl2004|) : esti 
mati ons for reactions involving S- a nd Si-bearing species 
from IWillacv fc Cherchnefll lll998[): combust ion mecha- 
nisms suchas^3^^Mec^JSmit^^^3^13) or the one 
of A. Konnov l|Konnovll2000|) : and from a revision of the 
last published data of rate constants for the type of re- 
actions mentioned above. When no data are available 
for the whole temperature range studied here (10-4000 
K) we have either extrapolated the expression or fixed 
the rate constant to its value at the nearest temperature 
for which it is known. When rate constants for reverse 
reactions were unknown or uncertain, detailed balance 
has been applied for calculating them from thermochem- 

1 The NIST Chemical Kinetics Database is available on the 
World Wide Web at http://kinetics.nist.gov/index.php 



TABLE 2 

Reactions relevant for oxygen chemistry 



Reaction A 




n 


7^ — 

c 


Important reactions in the lNNlOR, 


cnvelopc 




D1 U 1 U 1 tT . TT t TJ Q a w 1 fi- 
ll 1 H -f- 11 -j- ri2 — * it -j- 112 o.DXIU 


3.3 


-0.60 


u 


pn tt 1 it 1 TJ . TJ 1 T-T Q Qvin- 
ixZ a + ±1 + ri — > rl -+- n 0.0XIU 


33 


0.00 


U 


DO TT 1 TJ 1 TT , TT 1 TJ„ 1 w 1 f) — 

rto rl + ±1 + rle — > 11 2 H- tie O.IXIU 


33 


-0.13 


•~->t~\ 


T> A TT 1 TT TT 1 TT TT Q A -i r\ — 

H4 H2 ~~r -H-2 — > H + H ~r 11 2 0.4X1U 


09 


-0.24 


OZU4v5 


DC U IT , TT 1 TJ 1 TT Q ft w 1 O - 

ixD t±2 T~ tl — > 11 11 -j- rl o.DXIU 


09 


0.36 


0ZU4o 


TT ft TJ 1 TJ , JT 1 TJ 1 TJ „ C (Twin - 

ito H2 + rle — > 11 n -+- lie o.yxiu 


09 


0.23 


OZUUo 


rt( ol + OU — > blU + C l.oXlU 


09 


0.00 


34513 


R8 01O + (J — > 01 + CO 1.0x10 


09 


-0.23 


1291 


"D n /" 1 TT TJ 1 TJ 1 , , 1 r\ — 

K9 O + -H-2 —> OH + H 3.1x10 


10 


0.16 


11894 


rim r\TT 1 ■ c ■ f~~\ iTT 1 n . . 1 r\ — 

R10 OH + Si — > SiO + H 1.0x10 


10 


0.00 





T» 1 1 l~\ T T 1 TT T T f~\ iTT OO..ir\ — 

Rll OH + H2 — » H2O + H 2.2x10 


12 


1.43 


1751 


R12 biO + CO — * CO2 + 01 4.2x10 


13 


0.67 


32225 


T» 1 O /~1 i^"\ 1 T T TT /~\ 1 /*\ n , , 1 n- 

R13 OO2 + H2 — > H2O + OO 3.2x10 


7 


1.53 


56906 


K14 OH + OO —* OO2 + H 1.2X1U 


13 


0.95 


-73 a 


ni r /~\ 1 TJ /~\ TTiTT O 

H15 O + H2 — > OH + H 3.0X10 


13 


2.60 


3241 


R16 CO2 + Si — ► SiO + CO 2.7x10 


11 


0.00 


282 


Important reactions in the OUTER, 


envelope 




017 4- MWo i OTT -L- 1\TTT 1 9vin~ 

TX.l / V./ -f- 1N-TT2 " — ' V T IT -\- iNrl 1.ZX1U 


11 


0.00 


n 
u 


r 1 q 4- PRi — + Horn 4- h 1 4vin~ 

TT TO V./ ~p TT J * TT 2 \_s ~p TT 1.4: A TU 


10 


0.00 





R19 OH + S — * SO + H 6.6x10" 


11 


0.00 





R20 O + SH — * SO + H 1.2x10" 


10 


0.00 


-74 c 


R21 O + NH -4 NO + H 1.2x10" 


10 


0.00 





R22 CN + OH -* NCO + H 1.4x10" 


10 


0.00 





R23 S + CH 3 -> H 2 CS + H 1.4x10" 


10 


0.00 





R24 O + OH — > O2 + H 1.8x10" 


11 


0.00 


-175 b 


R25 O + CH2 — * HCO + H 5.0x10" 


11 


0.00 





R26 13 C + + 12 CO -> 12 C+ + 13 CO 3.4x10" 


10 


-0.50 


1.3 


R27 SiO + C + -> Si + + CO 5.4x10" 


10 


0.00 





R28 H^" + H2 — * + H 2.1x10" 


9 


0.00 





roq tj+ 1 nn TTr i rj+ _u T-Tr, 1 7vin~ 
nzy rig -\- * nuu ^ 112 x. 1 x iu 


9 


0.00 


u 


R^n 74(^0+ -u ttpn — * ro 4- Hr i T\ri4+ 1 1 Yin- 

rv.ou tt ^ yj -f- rToi n r \_j -\- tt\_^in tt o.ixiu 


9 


0.00 


u 


pqi TJ/~ir\-\- 1 „ jj„ c^c\ 1 , jj + 1 ;i vin- 
rlol HOVJ + U21I2 — > tU + L^2H^ 1.4X1U 


9 


0.00 


U 


ItoZ olo + — > olo + O Z.oXlU 


9 


0.00 


U 


dqq c;c+ 1 tj > CTj 1 c;+ 1 Q w 1 n~ 

rloo olo -f- n — > oil -+* ol 1 i.yxiu 


9 


0.00 


U 


R^4 SiS 4- S+ — * SiS+ 4- S ^ 9-xrin~ 

TTOT 010 O ' OlO O O.^AlU 


9 


0.00 





R35 SiO+ + CO — > CO2 + Si+ 7.9x10" 


10 


0.00 





R36 O + H -» OH + h^ 9.9x10" 


19 


-0.38 





R37 S + CO — > OCS + h^ 1.6x10" 


17 


-1.50 





R38 HCO+ + e- -» CO + H 1.1x10" 


7 


-1.00 





R39 H3CS+ + e~ — > H2CS + H 3.0x10" 


7 


-0.50 





R40 H 2 CO+ + e" -> HCO + H 1.0x10" 


7 


-0.50 





R41 SiO + hv -» Si + O 1.0x10" 


10 


0.00 


2.3 


R42 H 2 + hu -» OH + H 5.9x10" 


10 


0.00 


1.7 


R43 CH 4 + hv -> CH 3 + H 2.2x10" 


10 


0.00 


2.2 


R44 H 2 + CR -* H+ + e- 1.2x10" 


17 


0.00 





Reactions not considered in previous chemical models d 


R45 O + H 2 -> H 2 + hi/ 1.0x10" 


18 


0.00 





R46 HCO+ + H 2 -» H 3 CO+ + hi/ 5.0x10" 


IS 


0.00 






Note. — Rate constants for termolecular and bimolecular re- 
actions are given by k=Ax (T/300)" xexp(-C/T), for photodis- 
sociations (hu) are given by k=Axexp(-C A„) while for reactions 
with cosmic rays (CR) are given by k=A. Units are on 6 s"' for 
termolecular reactions, cm? s^ 1 for bimolecular reactions and s _1 
for photodissociations and reactions with cosmic rays. 

a k(T<80K)=k(T=80K) 

b k(T<167K)=k(T=167K) 

c k(T<298K)=k(T=298K) 

d See text for a discussion on these rate constants. 

ical properties of the species involved. Photodissocia- 
tion rates have been taken from the quoted databases 
when available, or assumed to be equal to those of sim- 
ilar molecules otherwise to ensure that all molecules are 
dissociated at some radius within the envelope. CO and 
H2 are known to be affected by self-shielding against in- 
terstellar UV photons due to their large abundance. CO 
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TABLE 3 

The first 5 reactions of the full chemical kinetics mechanism 



N 


Reaction 




A 


n 


c 


Ti 


T 2 


Trange 


Ref 


1 


H + H + H2 <=> H2 - 


f H2 


8.640E-33 


-0.60 


0.00 


50. 


5000. 


ALL__T 


NIST 


2 


H + H + He <=> H2 - 


1- He 


6.070E-33 


-0.13 


-39.00 


77. 


2000. 


ALL T 


NIST 


3 


H + H + H <=> H2 + 


H 


8.820E-33 


0.00 


0.00 


50. 


5000. 


ALL T 


NIST 


4 


CH + He <=> C + H 


4- He 


3.160E-10 


0.00 


33704.00 


2500. 


3800. 


LOW T 


NIST 


5 


CH + H2 <=> C + H 


+ H2 


8.830E-10 


0.00 


33704.00 


2500. 


3800. 


LOW T 


NIST 



Note. — The full version of tablc[3]is available in the electronic edition of the Astrophysical Journal. 



photodissso ciation rate as a funct ion of radial position is 
taken from iDotv fc Leund l)1998|) . who studied the CO 
case specifically for IRC+10216. For H2 we use the same 
result correcting for the different unattenuated photodis- 
sociation rate. 

The most important reactions for oxygen chemistry are 
shown in table El ■ Table contains the complete set of 
reactions. 

4. RESULTS AND DISCUSSION 
4.1. Oxygen chemistry in the inner envelope 

In this section we (i) first compare the results of chem- 
ical kinetics (assuming an expansion at constant velocity 
of 1 km s _1 ) with those of chemical equilibrium, (ii) sec- 
ondly we discuss the effects that SiO depletion from the 
gas phase could have on oxygen chemistry and (iii) finally 
we consider the chemistry with shocks. 

(i) The densities considered in the inner envelope af- 
fect, among other parameters, the H/H2 ratio. Figs.^ 3 , 
HI; and^i show the H2 and H abundances for models MH, 
MI and ML. The high densities of model MH make the 
chemical time scale to be lower than the dynamical time 
scale associated to an expansion at 1 km/s; therefore, 
reactions Rl to R6 

R1,R2,R3 H + H + M -> H 2 +M 
#4, i?5, R6 H 2 +M — > H + H + M 

are rapid enough to produce LTE abundances for H and 
H2. The steep decay in density of model MI makes that 
at ~2 R* reactions Rl to R6 become too slow, so that 
the H abundance is frozen at a value of 5xlCP 2 . Both 
models MH and MI allow hydrogen to enter the interme- 
diate shells and outer envelope mostly as H2. Model MI 
has a low photospheric density and therefore hydrogen 
exists in atomic form rather than molecular. For larger 
radii the decrease in temperature makes molecular the 
most stable form. However, chemical kinetics does not 
allow for a H^H2 transformation and hydrogen enters 
the outer envelope in atomic form. It is seen how both 
the density at the photosphere and the density radial 
profile affect the H/H2 ratio. There are observational 
constraints about this ratio that suggest that hydrogen 
is mainly molecular in the CSE (based on an upper limit 
for the H2 mass in the CSE and on an estimation of the 
H mass from observatio ns of the 21 cm line, see discus- 
sion in Glassaold 1996). This suggests a high density 
scenario, such as the one of models MH or MI, for the 
inner envelope. 

In what concerns oxygen chemistry, according to chem- 
ical equilibrium, for the densities and temperatures pre- 
vailing in the inner envelope, almost all the oxygen avail- 
able in a C-rich gas is locked in CO, with SiO being 2-5 
orders of magnitude less abundant. Any other O-bearing 



molecule has an LTE abundance always below ^10~ 10 . 
But some stellar radii away from the photosphere, where 
the temperature has decreased below ~700 K, the chemi- 
cal system experiences a significant change and molecules 
such as H2O and CO2 become very abundant in LTE (sec 
dotted lines in Figs. 0k, Eb and|2t ). This is related to 
the fact that the most stable form for both carbon and 
oxygen is the CO molecule in a high temperature regime, 
whereas for a low temperature regime carbon is p refer- 
ably locked in CH 4 and oxygen in H 2 l)Tsuil20 00'l. as is 
easily seen from the analysis of the equilibrium constant 
of this reaction: 

CO + 2,H 2 ^ CHi + H 2 (9) 

which can be fitted by the expression 

/ T \ 09 

K eq {T) = 1.58 x KT 50 ^) exp(24427/T) (10) 

for a temperature interval of 200-2000 K (K eq with units 
of (molecules/cm 3 ) -2 ). The formation of CH4 and H2O 
from CO and H2 is exothermic and involves less particles, 
so according to Le Chatelier's principle the CH4/H2O 
mixture dominates over the CO /H2 mixture at low tem- 
peratures and high densities. In our case, as temperature 
decreases so does density. Thus, a total CO— >H 2 trans- 
formation is only produced for T<400 K (r> 18 R*) in 
model MH (see Fig.EJi) and for larger radii in models MI 
and ML. Despite a total CO— >H20 transformation takes 
place beyond ^20 R* , the three models show a consider- 
ably increase in the LTE abundance of H2O beyond ^10 
R* . 

In order to have LTE abundances in the inner envelope, 
the chemical time scale T c hem has to be shorter than the 
dynamical time scale Td yn associated to the expansion for 
the gas being able to readapt its molecular abundances to 
the actual physical conditions (temperature and density) 
as it expands. The situation in the inner envelope of 
IRC+10216 is such that T c i iem <Td yn in the vicinity of the 
photosphere, but as temperature and density fall with 
radial distance, chemical reactions become slower and 
for a certain radio (that depends on the chemical species) 
i~chem>Tdyn- There, molecular abundances "freeze" and 
the LTE scenario disappears. 

The LTE abundance of CO remains constant from the 
photosphere to ^20 R*. Thus, whether chemical kinetics 
is rapid (LTE is reproduced) or slow ("freezing" effect) 
the CO abundance won't be modified. The situation 
is different for SiO (see Figs. EJi, Et> andEt). Its LTE 
abundance at ro is ^10 -8 while at ^3-5 R* it increases 
to ^10 -5 . For this abundance enhancement to occur as 
the gas expands, SiO has to take the oxygen from CO 
via reaction R7 

R7 Si + CO -> SiO + C 



(a) MH: High Density 
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(b) MI: Intermediate Density 



(c) ML: Low Density 
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Fig. 3. — Abundances, relative to total number of hydrogen nuclei, of some oxygen bearing molecules, —(a), (b) and (c): abundances 
given by LTE (dotted and dashed lines) and by chemical kinetics (solid lines) assuming a constant velocity expansion of 1 km/s and the 
density profiles MH— >(a), MI— >(b) and ML— »(c). —(d): abundances of CO and H2O for a CO,H2,CH4,H20 system with fc<j (see text) 
equal to 10 — 42 and 10 -46 cm 9 s . —(e) and (f): abundances given by chemical kinetics with the density profile MH when we consider that 
90% of SiO gas is deposited onto SiC grains and that the corresponding oxygen is released to the gas phase either as atomic oxygen— »(e) 
or as CO2— >(f). The nucleation rate of SiC grains J* (with units of nuclei cm -3 s -1 ) and its abundance xgic gr , relative to total number 
of hydrogen nuclei, are also plotted as dotted lines. See the electronic edition of the Journal for a color version of this figure 



which has a low rate and the above considerations about 
Tchem and Td yn apply. Only high densities and/or low 
expansion velocities will make SiO reach this abundance 
of ^1CP 5 . For the three density models considered and 
the expansion velocity of 1 km/s assumed, only the MH 
model makes SiO reach an abundance near its LTE value, 
whereas MI and ML lead to a non-LTE situation in which 
SiO abundance remains low. The abundance of atomic 
C has a great impact on SiO. In models MI and ML it 
remains well above its LTE value while in MH it follows 
the decreasing LTE profile being destroyed by reaction 
R9 

C + H 2 -> CH 



via reactions R10 and Rll 

R10 OH + Si -> SiO + H 
Rll OH + H 2 -> H 2 + H 

and (2) through CO2, directly via reactions R12+R13 



R12 
R13 



SiO- 

co 2 



CO 
-H 2 



co 2 - 

H 2 



Si 
CO 



R9 



H 



or involving OH which is converted into CO2 via reaction 
R14 

RU OH + CO -> C0 2 + H 
The following scheme shows these relations: 



mo 



Rll 



The H "freezing" effect of models MI and ML makes the 
reverse of reaction R9 to proceed at a higher rate that in 
a LTE situation making C overabundant which, through 
reaction R8 



R8 



SiO + C 



Si + CO 



destroys SiO making it underabundant with respect to 
a LTE situation. Again a high density scenario is re- 
quired for explaining th e SiO observational abun dance of 
8xl0~ 7 (estimated bv iKeadv fc Ridgwavl 1 199.il through 
IR ro-vibrational lines) that could be due either to this 
" freezing" effect or to depletion on SiC grains after reach- 
ing its LTE abundance of ~10 -5 . Neither of the models 
run made any other O-bearing molecule to have a signif- 
icant abundance. In this high temperature regime H2O 
is related to SiO by two routes: (1) competition for OH 



SiO OH — » H 2 

R12\ J.JH4 /R13 

co 2 

where the connection between the species is produced 
by the reactions shown above and also by their reverse 
processes, so we have a reversible system. Nevertheless, 
these mechanisms are not efficient in producing H2O 
from SiO. For H2O to take the oxygen from CO the 
chances are scarce because when the LTE abundance of 
H2O increases, the temperatures (<700 K) make that 
the possible chemical routes for a CO^H20 conversion 
will be too slow due to the large activation energy for 
breaking the CO bond. At this point, we would like to 
know if we miss some important reactions, rapid enough 
for allowing this CO^^O transformation in the condi- 
tions of the inner envelope of IRC+10216. We consider 
an scenario given by model MH (high density in order to 
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have a short r c ^ em ) and the chemical reaction given by 
equation 9. This reaction can be decomposed in several 
elementary reactions, for example 



CO + H 2 
OH + H 2 
CH + H 2 



fcs 



2 x 



CH 2 
H + H 



H 2 
H 2 
H 2 



H 2 
CH 2 
CH 3 
CH A 



OH 
H 
H 

-H 
-H 



(S) 



CO + 3 H 2 



H 2 + H 2 
CH 4 + H 2 



(G) 



In any mechanism composed of elementary reactions, 
the rate of the global reaction (G) will be given by the 
rate of the slowest elementary reaction (S). In this case 
the slowest elementary step will be the break of the CO 
bond, since it has the largest activation barrier. The 
number of H2O molecules produced or the number of 
CO molecules destroyed in the global reaction, per unit 
time and unit volume, is equal to the number of CO 
molecules destroyed in the slowest elementary reaction, 
per unit time and unit volume. We assign a rate con- 
stant ka for the global reaction G, calculate the reverse 
rate constant via detailed balance, and run the model 
for a CO,H2,CH4,H20 system in which only the global 
reaction G and its reverse operate. Then we vary kc 
until some significant fraction of CO is transformed into 
H 2 and CH 4 . Models for k G equal to 1CT 42 and 1(T 46 
cm 9 s _1 are shown in Fig. [2JI. It is seen that above a 
threshold value of ~1(T 44 cm 9 s" 1 for k G , the CO^H 2 
transformation is significant. The rate of destruction of 
CO molecules given by reaction G is kcxxco xni (xco 
is the CO abundance relative to H 2 and n is the density 
of H2) while the rate of destruction of CO molecules 
by reaction S is ksXxco xn2 ■ Since both rates have 
to be equal, we arrive to fcs^gxn 2 . Taking 10~ 44 for 
k G and a density of 10 10 cm~ 3 (roughly the value in 
the region 10-20 R* where the discussed process should 
take place) results in fcs=10 -24 cm 3 s _1 . Therefore, 
we need a bimolecular reaction similar to reaction S 
with a rate constant greater than 10 -24 cm 3 for 
temperatures lower than ~700 K. There is no data in 
the literature about a rate constant for a reaction such 
as S. Nevertheless, we can estimate an upper limit for 
ks if we assume a maximum temperature independent 
value of 10 -9 cm 3 s _1 for the reverse reaction and apply 
detailed balance. This results in activation energies of 
^80,000 K for ks, which for temperatures below 700 K 
imply rate constants well below 10 -24 cm 3 s _1 . So we 
don't expect a chemical route that is able to convert CO 
into H 2 in the inner envelope within the gas phase and 
during an adiabatic expansion. 

(ii) Such a CO^H 2 process seems too slow, but there 
is another O-bearing molecule which locks some impor- 
tant amount of oxygen, i.e. SiO, which is expected to 
condense onto SiC grains. How could H 2 abundance 
be affected by SiO condensation during the expansion?. 
When we consider such a process in our model, H 2 
abundance remains low. The reason is that, as explained 
above, H 2 is related to SiO in such a way that H 2 
becomes abundant when SiO does, but not when SiO 
depletes. Then, we consider in model MH (where SiO is 
formed with a significant abundance) a speculative pro- 



cess that can be summarized as follows: a SiO molecule 
condenses onto a SiC grain and the Si atom incorporates 
into the grain lattice while the O atom is released to the 
gas phase 

SiOgas + grain — > grain ■ ■ ■ Si + O gas 

We assume that 90% of the SiO formed in the gas phase 
is deposited onto SiC grains and that all that oxygen, 
formerly contained in SiO, is released to the gas phase as 
atomic oxygen (see Fig.Et). It is seen that molecules like 
H 2 and C0 2 increase their abundances. H 2 takes effi- 
ciently the oxygen when it is in atomic form via reactions 
R15+R11 



Rib 
Rll 



+ H 2 
OH + H 2 



OH + H 
H2O + H 



although some important fraction of the atomic oxygen 
returns to SiO from OH via reaction R10 



M0 



OH + Si 



SiO + H 



and through C0 2 via reactions R14+R16 

RU OH + CO -> C0 2 + H 



R16 



C0 2 + Si 



SiO + CO 



It is known in materials science that oxidation of solid 
SiC through a flow of 2 at temperatures of ^1000 K pro- 
duces solid Si0 2 and also CO and C0 2 molecules that dif- 
fuse trou£h_thelattice and are released to the gas phase 
(see lWang et al.l2001l and references therein). Therefore, 
some fraction of the oxygen entering solid SiC grains be- 
gin to build the Si0 2 lattice while the rest is released as 
CO and C0 2 . In the region of the CSE where SiC grain 
formation occurs, the temperatures are similar to those 
described above but the conditions are somewhat differ- 
ent since it is an oxygen-deficient environment. Never- 
theless, if SiO depletes on SiC grains, some fraction of 
the oxygen can be processed in the grain and released as 
CO and C0 2 . We have run our model assuming that the 
oxygen of the SiO depleted on SiC grains is released to 
the gas phase as C0 2 instead of atomic oxygen. Both 
H 2 and SiO compete for the oxygen contained in the 
C0 2 molecules through reactions R13 and R16 respec- 
tively 



R13 
R16 



co 2 
co 2 



H 2 
Si 



H 2 
SiO^ 



-CO 
CO 



but the high activation energy of R13 prevents water 
from reaching a high abundance. 

(iii) Now we will consider the effect that shocks have 
on oxygen chemistry. In Fig. 0] we plot the evolution of 
some abundances as the gas moves along the temperature 
and density profiles gi ven by model MH (see Fig. [2J , in 
order to compare with lWillacv fc Cherchnel l)1998jT who 
used high density values. We also plot the LTE abun- 
dances expected for the temperature and density profiles 
used, so we can see how the abundances predicted by 
chemical kinetics deviate from the LTE ones. Within 
the time interval associated to one shock plus relaxation, 
the immediately post-shocked gas has dissociated a cer- 
tain fraction of CO releasing atomic oxygen which has 
a large abundance during the first stages. This atomic 
oxygen is progressively returned to CO during the re- 
laxation time, but also to OH which has an apprecia- 
ble abundance during these first stages. Reaction Rll 
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Fig. 4. — Evolution of the abundances, relative to total num- 
ber of hydrogen nuclei, of some oxygen bearing molecules for 
a (a) 5 shocks history (N s f locka =5) and (b) 15 shocks history 
(N 3 f loc i eg =l5). The history starts at rg (i=0) and reaches r c at 
time NshocksXP- The LTE abundances for the physical conditions 
of the gas at each instant are shown as dotted lines. 

(OH+H2^H 2 0+H), which forms water from OH, oper- 
ates, but at these high temperatures the reverse reaction 
immediately destroys the H2O formed, returning oxygen 
into OH and some time later into CO. Indeed, shocks are 
a very good mechanism for dissociating CO and releasing 
atomic oxygen. This CO^O process occurs preferably 
during the first shocks due to the high temperatures, 
which also make the kinetics of the O^H20 process fast 
but, since water is not thcrmodynamically abundant at 
high temperatures in a C-rich gas (it is at low tempera- 
tures), other reactions are faster and return the atomic 
oxygen to CO. 

SiO is very affected by the high non-equilibrium chem- 
istry driven by shocks and ends this phase with a very 
low abundance, mainly due to reaction R8 , which was 
not included bv lWillacv fc Cherchnel j 1998ft . Note that 
either considering 5 or 15 shocks to reach r c from ro, the 
abundances do not appreciably change. It is shown that 
this mechanism decreases rather than enhances oxygen 
chemistry. 

Concerning the H/H2 ratio, shocks favor a large value 
since H 2 is dissociated in the post-shocked gas (due to 
the high temperatures reached) and depending on the 
density it can regenerate during the relaxation time via 
the three body reactions Rl, R2 and R3 (with high 
densities, see Fig. |5Ji) or freeze its abundance at some 
moment during the relaxation (with low densities, sec 
Figs. |2k and|5f). Again, with a shocks mechanism, a 



high density scenario is required to have hydrogen in 
molecular form. 

In summary, O-bearing species are not produced effi- 
ciently by gas phase chemistry in the inner envelope, nei- 
ther considering an expansion at constant velocity nor as- 
suming shocks. Nevertheless, the high water abundance 
obtained under the LTE assumption suggests that grains 
may be playing an important role in the recycling of oxy- 
gen from CO to H2O. 

4.2. Oxygen chemistry in the outer envelope 

The chemistry in the outer envelope is mainly driven 
by the photodissociation of molecules by the interstel- 
lar UV field followed by rapid neutral-neutral reactions. 
This mechanism forms species such as polyynes and 
cyanopolyynes from phot odissociation of p arent species 
such as C 2 H 2 and HCN ijMillar fc HerbsH 17994^ . Con- 
cerning oxygen, the photodissociation of CO occurs at 
the very outer edge of the envelope due to self-shielding 
(much further away than the photodissociation of C2H2 
and HCN) so atomic oxygen is abundant only at large 
radii in the CSE, where the density has decreased con- 
siderably and reactions are very slow. This fact makes 
oxygen chemistry to be much poor than carbon chem- 
istry, since oxygen keeps locked into CO along most of the 
envelope. The inclusion of the isotopo logue 13 CO, ~45 
times less abundant than parent CO ijCernicharo et all 
2000), which does not self-shield due to its lower abun- 
dance, can enhance O abundance at shorter radii. This 
allows for the formation of some O-bearing species with 
moderate abundances through rapid neutral-neutral re- 
actions and radiative associations. 

Ionic chemistry is triggered by cosmic rays ionization. 
Cosmic rays can penetrate deeper than UV photons in 
the CSE since they are not affected by opacity. Ions 
can react very rapidly via ion-molecule reactions but the 
chemistry they can induce is limited by the low abun- 
dance of ions in the CSE except at the very outer edge 
(where for example all the carbon is converted into C + ). 
In our model, ionic chemistry does not noticeably affect 
the abundances of the O-bearing species formed except 
for SO, whose abundance is greatly enhanced, and for 
HCO + , produced with a low abundance but enough for 
being detected due its large dipole moment (4.07 D). An- 
other important effect of ionic chemistry is the fraction- 
ation of CO through the exchange reaction R26 

R26 13 C++ 12 CO -> 12 C++ 13 CO 

which is known to be exothermic by 35 K. Below this 
temperature the reverse reaction does not operate and 
12 CO is effectively transformed into 13 CO. The net ef- 
fect is a delay in the destruction of 13 CO: its abundance 
remains nearly constant due to a steady state in which 
it is formed by the mentioned exchange reaction while it 
is destroyed by UV photodissociation. 

The abundances predicted in the outer envelope for 
CO, 13 CO, O, SiO and HCO+ (usual O-bearing species 
in a C-rich star) are shown in the left part of Fig. 03 . The 
second column of panels show the abundances predicted 
for H2O, OH and H2CO (the three "unexpected" oxy- 
gen species), while the rest of panels to the right show 
the abundances for some other O-bearing species. H2CS 
abundance distribution is also shown together with those 
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of C2H and CN, peaking at ~15" and ~20" respectively 
(see appendix B). We now describe individually the reac- 
tions that lead to their abundances, which are expressed 
relative to total number of hydrogen nuclei rather than 
relative to H2: 

- SiO is formed in the inner envelope and its abun- 
dance remains nearly constant throughout the outer en- 
velope until it is dissociated by UV photons (R41). Ionic 
chemistry destroys SiO at a shorter radius by reacting 
with C+ via R27 

R27 SiO + C+ -► Si+ + CO 

However, we point out that SiO could enter the outer en- 
velope with a much lower abundance due to condensation 
on grain surfaces. 

- HCO+ is produced by the same sequence of reactions 
that forms it in dark clouds. 

i?44 H 2 + Cosmic Ray — > H£ + e~ 
R28 H+ + H 2 -> H+ + H 

R29 H+ + CO -> HCO+ + H 2 

First, H2 is ionized by cosmic rays (R44), then Hj re- 
acts with H2 producing Hjj" (R28), which in turn re- 
acts rapidly with CO through a proton transfer reaction 
(R29). The origin of HCO + is different with respect to 
the rest of oxygen bearing species commented in this sec- 
tion because its formation does not require the photodis- 
sociation of CO but the presence of H^", the formation 
of which is related to cosmic rays. HCO + is abundant 
at shorter radii than the rest of O-bearing species with 
a peak abundance of 4.5xl0~ 10 at 3xl0 16 cm, and a 
corresponding column density of 1.1 xlO 12 cm~ 2 . The 
abundance profile sharply falls after its peak mainly due 
to dissociative recombination R38 (electron density be- 
comes important beyond 5x 10 16 cm) and to a less extent, 
due to proton transfer with HCN and C2H2 (reactions 
R30 and R31 respectively). 

i?38 HCO+ + e- -> CO + H 

R30 HCO+ + HCN -> CO + HCNH+ 

R31 HCO+ + C 2 H 2 -> CO + C 2 H+ 

H2O formation in C-rich environments, where 
atomic oxygen is available from the photodissociation of 
CO, is possible through the route of reaction R15 fol- 
lowed by reaction Rll. 

R15 + H 2 -> OH + H 
Rll OH + H2 -> H 2 + H 

This is the mechanism that form s water in the C-r ich pro- 
toplanetary nebula CRL618 fsee lCernicharol2004[) . where 
temperatures in the inner slowly expanding envelope are 
~300 K. However, in the outer envelope of IRC+10216 
such high temperatures are not expected and this route 
does not work, since both reactions R15 and Rll have 
activation energies that result in negligible rates. In our 
model, water forms through the radiative association be- 
tween atomic oxygen and molecular hydrogen (R45). 

i?45 + H 2 -> H 2 + hv 

To our knowledge, neither experimental nor theoreti- 
cal studies have made an estimation of this rate. We 
have assumed a temperature independent rate constant 
of 10~ 18 cm 3 s _1 , typical of radiative associations be- 
tween neutral species. Reaction R45 implies breaking 



the H2 bond so that atomic oxygen can insert in, which 
will result on some a ctivation energy . However, similar 
reactions, e.g. H9+C llLe Teuff et al.ll2000l) and H 2 +CH 
ijBrownsword et al]|l997|) take values 10- 17 -10- 18 cm 3 
s _1 . Water is produced over an extended region from 
4xl0 16 to 6xl0 17 cm with a peak abundance of 10~ 10 
and a column density of 4xl0 n cm~ 2 . These values 
are far below the abundance es timated from the ob ser- 
vations of SWAS (2-12xlQ- 7 . iMelnick et alJ 12001 or 
ODIN fl.2xlQ- 6 . [Hasegawa et al.l2006j) bv at least three 
orders of magnitude. We will discuss in §4.3 if water orig- 
inated in this way can or cannot explain the observations. 

- OH is produced from the photodissociation of water 
(R42) and from atomic oxygen through both the bimolec- 
ular reaction R17 

R17 + NH 2 -> OH + NH 

and the radiative association between O and H (R36). 
As the gas expands, reaction R36 becomes slower due 
to the decrease in density although it is compensated 
by the increase in atomic oxygen abundance. The re- 
sult is an increase of OH abundance only at the very 
outer envelope, shown in Fig. H as a peak at 6xl0 17 
cm. In fact, the O-bearing molecules formed by ra- 
diative associations show an extended distribution (re- 
actants such as atoms increase their abundances while 
expanding) whereas those formed by bimolecular reac- 
tions show a peak distribution (reactions only proceed 
in a shell where reactants are abundant before they are 
photodissociated). OH could also be produced in re- 
actions of O with several abundant hydrocarbons (e.g. 
C„H2, C„H, HC2n+iN) through H abstraction. How- 
ever, the reactions of atomic oxygen with closed-shell 
species such as C2H2, C4H2 and HCN have activation 
barriers whereas in reactions with radicals the H ab- 
straction channel seems to proceed slowly. For exam- 
ple, the reaction O+CH— >C +OH has been studied by 
iMurrell fc Rodriguez! {l986) who estimated a negligible 
reaction rate at low temperatures while the reaction of O 
with C2H produces mainly CO and CH. OH is predicted 
with a very extended distribution from 4xl0 16 to 10 18 
cm with a peak abundance of 7x 10 -10 at the very outer 
edge a nd a column density of 2xl0 n cm -2 . iFord et all 
( 2003) have derived an OH abundance of 4x 10~ 8 assum- 
ing that all OH is produced from the photodissociation 
of H2O, which is formed through comet evaporation in 
the inner envelope. The spatial distribution they obtain 
ranges from 3xl0 16 cm to 7xl0 16 cm. If OH follows 
the more extended radial distribution calculated by us, 
then the beam filling factors and physical conditions of 
th e emitting gas are very different from those assumed 
bv lFord et al.l l)2003t) and we could expect a considerable 
decrease of the OH estimated abundance. 

- H2CO is produced through reaction R18 

R18 + CH 3 -> H 2 CO + H 

the source of CH3 being the photodissociation of CH4 
(R43) while atomic oxygen is provided by CO and 13 CO 
photodissociation. The rate for R18 is well known for 
a temperature range of 259-2500 K in which it does 
not show a temperature dependence, suggesting a simi- 
lar value at low temperatures. The peak abundance is 
1.4xl0~ 9 , reached at 4x l0 16 cm. The abundance de- 
rived bv lFord et alJ l)2004[) from millimeter observations 
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Fig. 5. — Abundances, relative to total number of hydrogen nuclei, for some oxygen bearing molecules in the outer expanding envelope 
as given by chemical kinetics. C2H and CN abundances are shown in the bottom-right panel peaking at ~15" and ~20" respectively (see 
appendix B). Angular distance is given in the top axis for an assumed distance to IRC+10216 of 150 pc. 



of pure rotational transitions is a factor 5 higher than 
the one we get with the model. We discuss in §4.3 the 
compatibility of the abundance predicted by our model 
with the lines observed. 

- SO can be formed from OH+S (R19) and from 
O+SH (R20). 

R19 OH + S -> SO + H 
R20 O + SH -> SO + H 

For R19 we take the rate constant measured at 298 K. 
The rate of reaction R20 has been measured at 298 K and 
at high temperatures (1100-2000 K) and does not show 
temperature dependence, so that we assume the same 
value at low temperatures. In our model the abundance 
of OH is low in most of the CSE. On the other hand, SH 
abundance is relatively high due to a sequence of ion- 
molecule reactions that begin with the charge transfer of 
C+ to SiS (R32). The SiS+ formed this way reacts with 
H via reaction R33 and produces SH. 

i?32 SiS + C+ -> SiS+ + C 



- NO is formed by the bimolecular reaction R21 in- 
volving atomic oxygen and NH, which in turn comes from 
photodissociation of NH3. 



R21 



O + NH 



NO + H 



The rate for this reaction is a temperature independent 
estimation for the 250-3000 K range and is assumed to 
apply also at low temperatures. The predicted abun- 
dance is 2x 10~ 10 with a column density of 2x 10 11 cm~ 2 , 
although its low dipole moment (0.159 D) makes difficult 
its detection. 

- NCO forms by reaction R22 the rate of which is 
based on an experimental measurement at low pressure 
and 292 K and the assumption that formation of NCO 
is the main channel. 



R22 



CN + OH 



NCO + H 



R33 



SiS^ 



H 



SH + Si+ 



The peak abundance predicted for SO is 1.1 xlO -9 with 
a column density of 4.7xlO n cm' 2 . However, we point 
out that SO abundance is greatly influenced by H abun- 
dance, for which we hav e taken an upper limit given by 
iBowers fc Knannl l|1987f) . A process of H 2 formation on 
grains surfaces, not considered in this model, could re- 
duce the abundance of atomic H, thus decreasing the 
amount of SO in the outer envelope. 

- OCS forms by radiative association of S and CO 
(R37). The predicted peak abundance is 6xl0~ 12 with 
a quite low column density of 7.4xl0 9 cm~ 2 . 



Its predicted abundance is 10~ n . The relatively low 
dipole moment (0.6 D), the large partition function and 
the low predicted abundance will make very difficult its 
detection. 

H2CS abundance is also shown in Fig. Although 
it is not an O-bearing species, it is formed through the 
sulfur analogous reaction that forms H2CO (R23). 



R23 



S + CH 3 



H 2 CS + H 



The rate of this reaction has not been measured and 
therefore it is assumed equal to the oxygen analogous 
case. Atomic sulfur is necessary, instead of oxygen, and 
it is provided by the photodissociation of SiS and CS. 
Ionic chemistry slightly enhances H2CS abundance since 
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atomic S is released at short radii through ion-molecule 
reactions such as R34 

RM SiS + S+ -► SiS+ + S 

and because more chemical routes for thioformaldehyde 
formation appear, e.g. from HaCS" 1 " through dissociative 
recombination (R39). 

i?39 H 3 CS+ + e- -> H 2 CS + H 

The peak abundance is 10 -9 with a distribution similar 
to that of formaldehyde. The predicted abundance pro- 
file produces line shapes in reasonable agreement with 
observations (see below). 

- CO2 is formed through reaction R35, which involves 
SiO + that is formed with an abundance of ^lCP 11 from 
several ion-molecule reactions. 

R35 SiO+ + CO -» C0 2 + Si+ 

We point out that CO2 could be also a product of evap- 
oration of cometary ices and will have a large abundance 
in the cometary scenario. However, an ISO spectrum of 
IRC+10216 in the mid-infrared with a very large signal- 
to-noise ratio l)Cernicharo et aj]ll999l) only shows C2H2 
and HCN lines in the 13-16 //m range and no evidence 
forC0 2 . 

- O2 formation is strongly related to the presence of 
OH from which is formed through reaction R24. 

i?24 O + OH -> 2 + H 

Its spatial distribution follows that of OH with an abun- 
dance 2 orders of magnitude smaller. 

- HCO is produced with a quite low abundance 
(~10 -12 ) mainly from atomic oxygen by reaction R25. 

R25 + CH 2 -> HCO + H 

Dissociative recombination of H2CO + (R40) also con- 
tributes to HCO formation. 

4.3. Comparison with observations 

In order to compare with available observations we 
have performed Monte Carlo radiative transfer calcu- 
lations for H 2 CO, H ? Q HCO+, SO and H ? CS with 
the co de described in iGonzalez- Alfonso fc Cernicharo 
( 1993). The CSE has been simulated by a spherically dis- 
tributed expanding gas with the radial profiles for den- 
sity, temperature and abundance of each species taken 
from the chemical model. Line profiles and intensities 
for several transitions are then generated after convolu- 
tion of the radiative transfer results with the main beam 
of the select ed telescope . 

- H 2 CO. iFord et alJ <|2004fl have observed four pure 
rotational lines of ortho-H 2 CO: 2 1>2 — li 1 (140.840 GHz), 

21.1 - li,o (150.498 GHz), 3i, 3 - 2i, 2 (211.211 GHz) and 

31.2 — 2i ( i (225.698 GHz) with antenna temperatures of 
~30 mK, ~20 mK, ~20 mK and ^20 mK respectively; 
one line of para-H 2 CO: 3 0:3 - 2 , 2 (218.222 GHz) with 
TJJ~20 mK; and failed to detect the low excitation tran- 
sition of para-H 2 CO lo,i — 0o,o (72.838 GHz) at a sen- 
sitivity of ~10 mK. They also conducted 17" offset ob- 
servations for the 2i.i — l^o line with an average line 
intensity reduced by a factor 3-4, which nevertheless in- 
dicates an extended H 2 CO distribution. We have used 
the abundance profile obtained from our chemical model 



to predict the expected line intensities in order to see if 
the formaldehyde formed through these chemical routes 
can account for the observed lines. The collisional co- 
efficients, need ed for solvin g the statistical equilibrium, 
are taken from lGreenl 1^991) (corrected for H 2 as collider 
instead of He) and the normal ortho-to-para ratio of 3:1 
has been assumed. 

The resulting lines are less intense than the observed 
ones by a factor 4-10 depending on the line. Nevertheless, 
one has to keep in mind that the abundances resulting 
from such a simple chemical model are affected by un- 
certainties due to missing reactions, uncertain rates and 
also errors in physical parameters (kinetic temperature, 
mass loss rate, etc). Formaldehyde is formed from O and 
CH3 by reaction R18 

#18 + CH 3 -> H2CO + H 

so the abundances of the two reactants together with 
the rate of R18 directly affect the abundance that H 2 CO 
may reach. The rate constant for reaction R18 is known 
down to 259 K, but not for the much lower temperatures 
that are present in the outer envelope of IRC+10216. We 
could expect a higher rate constant at very low tempera- 
tures, a behavior experimentally observed in many other 
radical-molecule r eactions without activation energy (see 
Smit h~et al J 120041 and references therein). Another im- 
portant parameter is the branching ratio for CH3 pro- 
duction in methane photodissociation by t he interstellar 
UV field, for which ILe Teuff et all (1200?!) give a 15%. 
For the solar UV field values up to 44% are commonly 
used (see the chem ical model of Saturn's atmosphere by 
iMoses et ai]l2000ft . Although the energy distribution of 
the interstellar and solar UV field are somewhat differ- 
ent, we could expect larger uncertainties in the former 
one. When we increase R18 rate constant by a factor 
of 3 and the CH3 branching ratio by a factor of 2, the 
resulting H2CO abundance (see solid line in right-down 
panel of Fig.El) increases by a factor of ~6 with respect 
to the former model (that of Fig. O ■ A comparison of 
the spectral lines resultin g from such abund ance profile 
with the lines observed bv lFord et alJ l)2004t) give us sev- 
eral conclusions: (1) the line intensities agree quite well 
for the ortho-H2CO lines whereas for the para-H2CO the 
agreement is worse: the 72.838 GHz line is predicted a 
factor ^3 more intense than the observed upper limit 
and the 218.222 GHz line is predicted a factor 2-3 less 
intense than the observation. Uncertainties in the colli- 
sional coefficients, ortho-to-para ratio or in the tempera- 
ture and density profiles used could explain the discrep- 
ancies between model and observations. (2) The ratio of 
the integrated intensity between the 17" offset position 
to the centered one in the 150.498 GHz line is 0.87, in 
contrast with the value of 0.4 obtained from the observa- 
tions. Sources for the disagreement can be errors in the 
pointing, the uncertainty in the distance to the star or 
the depart from spherical symmetry in H2CO emission. 
(3) One last aspect to point out is the different shapes of 
the observed and predicted 1.3 mm lines of ortho-H2CO 
(at 211.211 and 225.698 GHz). The predicted lines have 
pronounced double-peaked shapes with almost no emis- 
sion at zero velocity, as corresponds to a shell with a 
certain size (diameter of ~40" in the model) observed 
with a smaller telescope beam (~11" for IRAM-30m at 
those frequencies) . On the other hand the observed lines 
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be compared with the observations reported by Ford ct al. (2004). IRAM-30m beam size has been considered and the lines have been 
smoothed to a resolution of 1 MHz. The H2CO abundances, relative to total number of H nuclei, are plotted in the bottom-right panel. 
-Solid lines: H2CO abundance resulting from the chemical model with two parameters slightly changed: R18 rate constant increased in a 
factor 3 and the CH3 branching ratio in methane photodissociation increased from 15% to 30%. -Faint dotted lines: H2CO produced from 
photodissociation of an unknown pa rent species , whic h would be released by sublimation of the mantles of comets, assumed to be present 
in the inner envelope. Best model in lFord et al J 120041) . -Hard dotted lines: H2CO abundance resulting from the chemical model including 
reaction R46 in order to increase formaldehyde abundance in the inner regions. 



show substantial emission at zero velocity and a less pro- 
nounced double-peaked profile. Neither an error of up to 
4" in the pointing observations nor a larger distance to 
the star of up to 200 pc (more will mak e IRC+10216 to be 
overl uminous for its type according to lCrosas &: Mentenl 
Il997[) will appreciab ly change the predicted line profiles. 

iFord et alJ l)2004D interpret the presence of formalde- 
hyde as a direct product from the photodissociation of 
an unknown parent molecule released in the inner en- 
velope by cometary sublimation. Such interpretation is 
supported by the observation of extended emission of 
forma ld ehyde towa r d the comet P/Halley (Meie r~et alJ 
Il993|) . IFord et all l)2004|) constructed a simple model 
with the following parameters: initial abundance of the 
parent species relative to H2 of 8xl0 -8 , unattenuated 
photodissociation rate of the parent species producing 
H 2 CO of A p =1.6x 10~ 10 s _1 and unattenuated pho- 
todissocation rate of formaldehyde Ah 2 co=7-81x 10~ 10 
s -1 , which reproduced the observed line intensity of the 
150.498 GHz line and the ratio (/ T* A dv off )/{$ T* A dv cent ) 
of 0.4 for this line. We have run a model with the same 
parameters to obtain the abundance distribution and the 



expected line profiles (light dotted lines in Fig. [fJJl . It is 
seen that the abundance distribution is more extended 
than the previous one although the 1.3 mm line profiles 
of ortho-H2CO still show an emission at zero velocity too 
low compared with the observations. The abundance dis- 
tribution of H2CO with a cometary origin is somewhat 
similar to the one we get with the chemical model, be- 
cause in both cases formaldehyde is produced in the re- 
gion where interstellar UV photons penetrate. 

From the above considerations it follows that an in- 
crease of the formaldehyde abundance in the inner re- 
gions of the CSE is needed for producing some emission 
at zero velocity and matching the observations. Forma- 
tion of H2CO in inner regions can be explained by hydro- 
genation of CO on grain surfaces followed by photodes- 
orption (species such as CH4, NH 3 and H 2 S are thought 
to be formed by similar processes). Since this possibility 
is very complex to model, we have investigated if some 
missing gas phase reactions can enhance formaldehyde 
abundance at short radii. At this point it is worth to note 
that HCO + has an abundance distribution peaking closer 
to the star than the one of H2CO (see Fig 01 and that 
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produces double-peaked spectral lines but with substan- 
tial emission at zero velocity (see Fig. 03 ) . As discussed 
for HCO + formation, ionic chemistry is initiated by cos- 
mic rays at shorter distances than photochemistry does 
since the CSE is not opaque to cosmic rays. Formalde- 
hyde can also be formed from H2CO + (through charge 
transfer with several species) and H 3 CO + (through dis- 
sociative recombination) but in our model these routes 
are not important since both H 2 CO + and H 3 CO + are 
not abundant enough. Looking at exothermicity we can 
examine possible missing reactions that can contribute 
to the formation of these two species (heats of forma- 
tion have been taken fro m NIST Chemistry Webbook 2 
and lLe Teuff et aLll2000ft . One possibility is the reaction 
+CO for which the only exothermic channels are the 
well known proton transfer giving HCO + and the radia- 
tive association which produces H 3 CO + . Another possi- 
bility is the reaction HCO + +H2 but all the investigated 
channels are highly endothermic due to the high stability 
of both reactants and only the radiative association pro- 
ducing H 3 CO + is exothermic by 123 kJ/mol. We find 
that including the reaction HCO++H 2 -> H 3 CO++Iw 
with a rate constant of 5xl0~ 15 cm 3 s^ 1 (R46), the 
formaldehyde abundance is enhanced at the inner regions 
producing rotational lines in agreement with the obser- 
vations, i.e. with significant emission at zero velocity 
(strong dotted lines in Fig. |5J) . All the lines except the 
one at 218.222 GHz of para-t^CO (predicted a factor 
2 less intense than observed) agree very well with the 
observations in intensity and profile. The required rate 
const ant is typical of other ion-neutral radiative associa- 
tions ( Le Teu ff et allkOOOl) . 

Finally, is it wor th to note that t hree of the four 
unidentified lines in iFord et aD l)2004|) can be assigned 
as follows: U150 i s the 7i. 7-61.fi transition of S i 13 CC at 
150385.281 MHz l|Cernicharo et al.lll99ll 1200(1 : U218b 
is the iV=23-22 J =45/2-43/2 transition o f CC 13 CCH 
at 218103.281 MHz <|Cernicharo et al]l200(ih : and U218a 
is one component of the 2 Hi/2 J=45/2-43/2 transition 
in the y 7 =l vibrati o nal sta te of C 4 H at 218 2 87.460 
MHz (iGuelin et all Jl987bt lYamamoto et alJ Il987t 
iCernicharo eTal. 2000). 



- H2O. The evidence for water presence in IRC+10216 
comes from the detection of the ortho transition li,o- 
l i at 556.936 GHz with the telescopes SWAS (T^-20 
mK. iMelnick et all l200H) and ODIN (T M b~50 mK, 
iHaseeawa et al.ll2006|) .' In order to model the radiative 
transfer of water throughout the CSE we consider the 
lowest lying 8 rotational levels in the vibrational ground 
state level (fi,^,^)— (0,0,0), in the first excited bend- 
ing mode (0,1,0) and in the first excited asymmetric 
stretching mode (0,0,1) . The ortho-H20-He collisional 
rates from lGreenl fl993). corrected for H 2 as collider, are 
adopted for all the transitions within the same vibra- 
tional level, whereas for ro-vibrational transitions colli- 
sions are not considered important. 

Excitation to the ^2 = 1 state occurs by absorption of 
A~6 /im photons whereas excitation to the ^ 3 =1 state 
corresponds to A~3 /im. Photons at these two wave- 
lengths are abundant within the CSE due to the pres- 

2 The NIST Chemistry Webbook is available on the World Wide 
Web at http://webbook.nist.gov/chemistry/ 



ence of the central star and dust, so these excited vi- 
brational levels can be easily populated. Furthermore, 
since our chemical model predicts water to be abundant 
in the outer CSE, where densities range from 10 5 to 10 2 
cm~ 3 (n cr i t ^10 9 cm~ 3 for the 556.936 GHz transition), 
collisional excitation is not effective to populate the l^o 
state and the dominant excitation mechanism is radia- 
tive pumping to the 1^2 =1 and ^ 3 =1 vibrational excited 
states followed by radiative decay to several rotational 
levels of the vibrational ground state, included the l^o 
state. 

In our chemical model water is formed through the 
radiative association R45 



i?45 



+ H 2 H 2 + hi 



with an assumed rate constant of 10 18 cm 3 s 



As 



noted in §4.2, the abundance reached with such a rate 
constant is not enough to explain the observed 556.936 
GHz line. Since the rate constant for this reaction is not 
known and water abundance approximately scales with 
this rate, we have fitted R45 rate to get a 556.936 GHz 
line intensity in agreement with the observations. We 
find that it is necessary to increase the rate constant up 
to ^10~ 15 cm 3 s -1 . The peak abundance obtained then 
is ^10~ 7 as shown in Fig. [7| . With such an abundance 
profile, we have performed Monte Carlo radiative trans- 
fer calculations. We show in Fig. the lines predicted 
for ODIN in different scenarios in which the fluxes at 
A=3 and 6 /im vary. On the left side the central star 
has been considered as a blackbody at T—T^ with and 
without dust throughout the CSE. We consider carbona- 
ceous grains with sizes betwe en 0.01 and 1 /tm and a 
gas-to-dust mass ratio of 500 (jKnapd 1 1985(1 . It is seen 
that the presence of dust enhances the 557 GHZ line in- 
tensity since grains absorb the photons from the star (its 
maximum is at A=2 /tm) and reemits at longer wave- 
lengths, thus enhancing the flux at A=6/im and exciting 
the bending mode of H2O. In f act, the IR spectrum o f 
IRC+10216 as observed by ISO (jCernicharo et al.lll999]l 
peaks at 10 /tm due to the presence of dust and approxi- 
mately corresponds to a blackbody of radius 20 R* with 
a temperature of ^500 K. We can therefore simulate the 
IR flux that water in the outer envelope would receive by 
putting such a blackbody instead of the star. In Fig.0we 
see how the 556.936 GHz line intensity increases when the 
flux at A 6 /im is enhanced by considering a blackbody 
with increasing temperature. These radiative transfer 
models show the importance of the IR pumping of the 
excited vibrational levels of H2O for the intensities of 
pure rotational lines in the ground vibrational state. 

In summary, assuming a rate constant of 10~ 15 or a 
few 10~ 16 cm 3 s _1 for R45 we can reproduce the water 
line observed by ODIN. Now the question is how reason- 
able is such a high rate constant for a radiative associ- 
ation between neutrals. This reaction is spin-forbidden 
according to Wigner-Witmer rules since the potential en- 
ergy surfaces (PES) of reactants and products in their 
ground electronic states, H 2 ( 1 S)+0( 3 P)^H 2 0( 1 A)+lw, 
are not connected adiabatically. For the reaction to pro- 
ceed, an intersystem crossing between PES with different 
total electronic spin has to occur during the reaction. A 
spin-forbidden reaction is likely to have a small rate con- 
stant, however some other reactions in which reactants 
and products do not adiabatically correlate have been 
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Fig. 7. Left: Water vapor abundance, relative to total number of H nuclei, when assuming a rate constant for R45 of 10 — 15 cm 3 s — 1 . 

-Right: calculated line intensities and profiles for the li,o-lo,i transition of ortho-H20 for different IR fluxes. ODIN beam size (~2') has 
been considered and the lines have been smoothed to a resolution of 1.3 km/s. 



found to be very rapid. For example the bimolccular re- 
action: C( 3 P)+C 2 H 2 ( 1 S)^C3( 1 I])-(-H2( 1 S) may have a 
rate const ant of several 10~ 1 cm 3 s _1 at very low tem- 
peratures ijClarv et alJ l20021 while the radiative associ- 
ation HS+( 3 S)+H 2 ( 1 S)^H 3 S+( 1 A )+h£ has a rate co n- 
stant of 7x10" 16 cm 3 s" 1 at 80 K l|Herbst et al] ll9891. 

A high rate constant for reaction R45 would also 
have consequences for the chemistry in dark clouds. If 
k.R45=10~ 15 cm 3 s _1 , and with a typical density of 10 4 
cm~ 3 , the chemical time scale to produce water from 
atomic oxygen would be as short as 3x 10 3 years. Chemi- 
cal models without this radiative association predict H2O 
abundances somewhat higher t han the observational up - 
per limits imposed by SWAS ijRoberts fc Herbstl [20021. 
Nevertheless, depletion on grains at late times would re- 
duce the water vapor abundance so that models could 
match the observations. In diffuse clouds, with lower vi- 
sual extinctions (A v ~l) and densities (n#~100 cm -3 ), 
reaction R45 would also produce water in time scales 
of 3xl0 5 years, with an abundance relative to H2 of 
^3xl0 -7 in steady state 3 . This value is 1-2 orders of 
magnitude hig her than the scarce observational estima - 
tions available (Soaans et al.lll998t iNeufeld et al] [20021. 
Nevertheless, in diffuse clouds a significant fraction of hy- 
drogen is not in molecular form, therefore the rate of wa- 
ter formation through reaction R45 would be lower than 
our estimation. Furthermore, since the observations of 
interstellar water ice bands do not discriminate between 
different types of clou ds, dense or diffuse, present along, 
the li ne of sight (e.g. iSchutte et alJll998l: iMoneti et all 
12001(1 , it is highly uncertain whether the amount of water 
as solid ice is important or not. 

There are other clues, apart from the 557 GHz line, 
that will help in understanding the water origin in this 
late-type C-rich star: (1) If water is formed in the outer 
envelope it will have a spherical shell-like distribution 

3 The H2O abundance reaches the steady state when the forma- 
tion and destruction rates become equal. If we assume that reac- 
tions R45 and R42 are the main formation and destruction routes, 
then n(H 2 0)/n(H 2 ) = [k R45 x Q n H2 ]/[G 5.9xl0- 10 exp(-7.1 
j4„)], where xo is the O abundance relative to H2 and Go is the 
UV field relative to the interstellar standard one (see appendix B). 
If we take A v =l, Go=l, xo=6xl0 -4 , fci}45=10 -15 cm 3 s _1 and 
«.(H 2 )=n£f /2=100/2 cm~ 3 (assuming hydrogen is mainly molecu- 
lar), we get n(H 2 0)/n(H 2 ) = 3xlCT 7 . 



while Fischer- Tropsch and cometary hypotheses will 
produce a solid sphere-like distribution since both 
require water to be formed much closer to the star. 
Herschel Space Observatory will be able to distinguish 
between these two possibilities since its beam size will 
be much narrower than those of SWAS and ODIN. 
(2) The absence of pure rotational and ro-vibrational 
lin es of water in the ISO far-IR spectrum analyzed 
by iCernicharo et al.1 l)199fiH discards the presence of a 
significant a mount of water in th e inner envelope. (3) 
As noted bv iRodeers &: Charnlevl ((2002), the cometary 
origin implies the presence of HDO with an abundance 
^0.06 % relative to H2O. HDO detection would demon- 
strate that water is not formed from material ejected by 
the star and the cometary hypothesis would be the only 
one reasonable. 

- HCO+. Its presence in IR C+10216 has been de- 
bated for years. Chemical models l|Glassgold et all 1986(1 
predicted it with an abundance high enough to allow its 
detection with available millimeter or submillimeter tele- 
scopes. On the other hand, observations did not clearly 
probe its existence. The J=l-0 line (with an intensity 
of T^=20 mK with the IRAM 30m telescope, see Fig.EJl 
often has been observed as a weak feature at the noise 
level (|Lucas fc Guelini ri9901. The J=4-3 line was ob- 
served as a weak feature of T* ~25 mK with the JCMT 
telescope bv lAverv et al.l (|1994|1 . However, they failed to 
detect the J=3-2 line at a noise level of ~10 mK, which 
was predicted with Tj^45 mK on the basis of a rota- 
tional temperature of 17 K obtained from the intensities 
of the J=l-0 and J=4-3 lines. These considerations led 
these authors to conclude that the feature at the J=4-3 
frequency was not arising from HCO + and they calcu- 
lated an upper limit for its column density of 1.4xlO n 
cm~ 2 . Our Monte Carlo radiative transfer calculations 
show that the excitation temperatures of each line are 
significantly different, thus a prediction of the expected 
intensities of other lines with a unique T rot should be 
taken with care. To calculate the line profiles of HCO + 
we have decreased by a factor of 2 the abundance pro- 
file obtained with the chemical model in order to match 
the observed J=T0 line. HCO + is formed through the 
sequence of reactions R44+R28+R29 and therefore its 
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Fig. 8. — Calculated line intensities and profiles for pure rotational transitions of HCO + , SO, ortho-H2CS and para-H2CS. All lines have 
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J=l-0 line observed. The model is plotted as dotted lines and the A 3 mm observations from|Ccrnicharo ct al. (2006) as solid lines. 



abundance depends basically on the rates of these reac- 
tions as well as on the destruction rates. The most un- 
certain of all of them is the ionization by cosmic rays of 
H2 (R44), which depends on the cosmic rays field in the 
surroundings of IRC+10216. The abundance of HCO + 
approximately scales with R44 rate. Therefore, a de- 
crease by a factor of 2 in the cosmic rays field would 
explain the decrease in abundance needed to match the 
J=l-0 observed line. 

The predictions for the first four pure rotational 
transitions are shown in Fig. |H1 It is seen that the 
excitation conditions are such that the J=3-2 line is 
expected with a slightly higher intensity than the J=4-3 
line, both being ~10 mK in T* A for the JCMT telescope. 
This sugges t s tha t the 25 mK feature observed by 
lAverv et alJ (^994) at the frequency of J=4-3 is indeed 
not arising from HCO + because the J=3-2 was not 
detected above a 10 mK noise level. Future long 
time-integrated observations of these two lines will 
definitively establish the excitations conditions of HCO + 
in IRC+10216. 

- SO. The abundance reached by SO in the chemical 
model is enough to produce some lines with intensities 



of a few mK. Fig. |H1 shows the Nj=2%-\-2 and 34-23 
expected line intensities together with a feature at 
99.300 GHz o f a spectrum, obtained with the IRAM 
30m telescope l|Cernicharo et al.ll2006|) , which could cor- 
respond to the 23-I2 transition. A higher signal-to-noise 
ratio for this spectrum is needed to claim detection. 
Long time-integrated observations at the frequencies of 
other SO transitions (e.g. the 34-23) would decide about 
the existence of SO in IRC+10216 formed trough the 
ion-molecule chemical sequence of R32+R33 followed by 
the neutral-neutral reaction R20. 

- H2CS. The chemical model predicts that this 
molecule is produced with a moderate abundance, both 
through neutral chemistry (from S+CH3, R23) and 
through ionic chemistry (from H 3 CS + , R39). The abun- 
dance reached is within the same order of magnitude than 
for H2CO, but less intense rotational lines are expected 
due to its smaller dipole moment (1.649 D versus 2.331 
D of formaldehyde). We have also assumed an ortho-to- 
para ratio of 3:1. Fig.IEjshows the expected line profiles of 
ortho and para H2CS, together with some features from 
a A 3 mm spectrum take n with the IRAM 30m telescope 
ijCernicharo et al1 l2006). The features have intensities 
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somewhat higher than the predicted ones, but with the 
right order of magnitude. 

5. CONCLUSIONS 

LTE chemical models predict that O-bearing species 
such as water are very abundant in the inner envelope 
of IRC+10216. However, a more realistic non-LTE ap- 
proach, based on chemical kinetics, indicates that the 
transformation of CO into H2O within the gas phase is 
not efficient in the inner layers, because of the high ener- 
gies required to break the CO bond. On the other hand, 
the increase of SiO abundance with radius, predicted in 
LTE, is possible because the reaction Si+CO^SiO+C is 
competitive in a high density scenario for the inner enve- 
lope. An alternative mechanism for water production out 
of the gas phase could be related to grain surfaces, which 
can act as a catalyst reducing the activation energy for 
either a CO— >H 2 Q process (due to Fischer- Tropsch catal- 
ysis according to lWillac"vll2f)04ft or a SiO^H 2 process 
(on SiC grains). 

Concerning oxygen chemistry in the outer envelope, 
the release of atomic oxygen due to CO photodissocia- 
tion allows for the formation of some O-bearing species. 
The abundance predicted for H2CO is a factor ~5 
lower than the observational estimation of iFord et aP 
(2004). Wether formaldehyde is formed by the reaction 
0+CH 3 ^H 2 CO+H or it is produced from photodissoci- 
ation of an unk nown parent specie s released from comets, 
as suggested bv lFord et alJ 1)2004(1 . the expected shape of 
the A 1.3 mm rotational lines disagrees with the observa- 
tions. An alternative source of H2CO is suggested to be 
the radiative association HCO + +H2^H3CO + +1iza The 
possibility of forming water in the outer envelope de- 
pends on the rate constant of the radiative association 
between atomic oxygen and molecular hydrogen. Non- 
local radiative transfer models show that a rate constant 
as high as 10 -15 cm 3 s _1 is needed to reproduce the 
556.936 GHz line profile observed by SWAS and ODIN 



telescopes. Quantum chemical calculations of the rate 
constants of the two radiative associations suggested in 
this paper, specially the H2+O reaction, could support 
or discard them as important reactions in astrochemistry. 

Other oxygen bearing species such as SO could ex- 
ist with abundances and excitation conditions which 
would produce rotational lines near the detection limit 
of IRAM-30m telescope. HCO + is observed with an in- 
tensity roughly in agreement with that predicted from 
chemical and radiative transfer models. Thioformalde- 
hyde formation is also predicted through a chemical route 
analogous to that of formaldehyde. This is supported by 
the agreement of the line profiles obtained from radiative 
transfer models with observations at A 3 mm. 

Although the evaporation of cometary ices could be 
the source of water in this late-stage carbon star, other 
phenomena such as catalysis on grains or the production 
in the outer envelope through the radiative association 
of H2+O could also explain the formation of this O- 
bearing species in the expanding carbon-rich envelope. 
Future observations with the Herschel Space Observatory 
will permit to distinguish the most plausible chemical 
processes leading to the formation of water vapor in a 
carbon- rich environment. 
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APPENDIX 

NUCLEATION AND GROWTH OF SIC GRAINS 

Here we describe the formation of SiC grains in the inner envelope of IRC+10216 as a two-step process: (1) formation 
of condensation nuclei plus (2) growth of grains from these nuclei by accretion of gas phase species. 

(1) The mathematical description of the nucleation process is mainly taken from iGail et al.l J1984). We assume 
that the expanding gas has at all moments a population of clusters (SiC) at, composed of different number N of SiC 
monomers, which in thermal equilibrium is given by 

n { siC) N = n SlC x exp( - AG(N)/kT gr ) (Al) 

where n means numerical density, T gr is the grain temperature and AG(N) is the free energy of formation of a cluster 
of size N from the monomers. The temperature of SiC grains T gr is not equal to the kinetic temper ature of the ga s 
Tfc because the former is affected by an inverse greenhouse effect which makes T gr to be less than Tk l|McCabel fl982'). 
To calculate T gr we follow the treatment of this author. The magnitude AG(N) can be expressed as 

AG(N) = k6 N (T gr )(N - I) 2 / 3 - kT gr (N - l)lnS (A2) 

where S is the supersatu ration ratio (ratio of the SiC vapor press ure to its sat uration pressure) that is calculated 
following iMEHiEi CHI, and the magnitude 0N(T gr ) (described in lDraine1ll979D is related to the surface tension of 
solid SiC for large N while for low N is treated as a free parameter since it is unknown. 

Once we have establish a population of clusters (SiC) n given by thermal equlibrium, nucleation theory says that 
there exist one critical size 2V*. Cluster with sizes N>N* will continue accreting growth species (here SiC molecules) 
resulting in grains of growing size, while clusters with sizes N<N* will revert to the monomers. Thus, there exists a 
bottleneck at A* in the size spectrum. The critical size A* is found by maximization of AG(N) with respect to N. 
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The steady state rate Jf for the transformation '(SiC)jv„— ► nucleus' is given by 

I N< 

An, 



Ei kT k \l/2 1 _ . 2 * r 2/3 



Xn (SiC) N , x ^ 



(A3) 



The subscript i refers to the growth species: is the numerical density, rrii the mass, a.i the sticking coefficient, and 
the sum of all contributions (3 has the meaning of a flux of particles. The subscript * refers to critical clusters: An, 
is the surface ar ea and ao is the radius of a monomer in the SiC lattice. Z is the Zeldovich factor (see Section 2 of 
iGail et al"]ll984l for a detailed explanation) . 

Once is computed for the conditions prevailing at each radius, the temporal evolution of the nucleation rate </* 
is obtained from just considering a time lag r* ( see IGail et al] IT9841 



d f J, \ I [ J*- Jf ] 



dtlpA 



(A4) 



(2) The growth of grains occur by addition of species that incorporate into the grain, resulting in a temporal evolution 
for the grain radius a gr of 

' kTk \ 1/2 (A 2 \ "II. 
-27TTOi 



da 



dt 



2^ a i X ) X ( 47r V) X X 



3A„ 



2/3 



(A5) 



where the sum extends to all gaseous species i that produces growth of grains (here only SiC molecules are considered). 
The parameters a, and m, have the same meaning as in equation A3, ao j is the radius of the species i in the lattice 
and Ni is the number of monomers of type i forming part of a grain. 

The unknown parameters, such as the sticking coefficients and the form of the function On for low N are chosen for 
obtaining an abundanc e and size of SiC grains in agreement with the literature: abundance relative to H2 of ^10~ 14 
and radius of ~0.1 /xm orenz- Martins fc Lefevrdll993() . 

VISUAL EXTINCTION RADIAL PROFILE 

The parameters Ai and Ci in equation 8 are usually given in the literature for plane-parallel geometry. Therefore, 
the evaluation of the photodissociation r ate Tj of a species i through equation 8 is only valid for a species situated at 
point P in a plane-parallel layer (see Fisr. lBll ). where the visual extinction measured along the direction perpendicular 
to the infinite plane is A^ p . But we want to know Ti at the different points r in spherical geometry. The strategy, 
then, is to find expressions for the UV field (AnJuv) at points P and r in plane-parallel and spherical geometries 
respectively. Equaling both expressions we will determine the associated A^ p , that we will insert in equation 8. 
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Fig. Bl. — A v ra dial profiles obtained from the plane-parallel to spherical correction (Ay P ) and from the standard relation A v ccNh of 
Bohl in et alJ 119781) . The inserts show a two dimensions scheme of plane-parallel and spherical geometries. 



The UV field at point P in a constant density plane-parallel layer is given by 



'Tin 00 \ App 1 



4irJ uv (P) = 2vrG / sin6exp{ - T 1000 (9)}d6 , t 100 o(#) = f — ) T 



500 ' 1.086 cosO 



(Bl) 
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where Go is the wavelength integrated specific intensity of the interstellar standard UV field l)Drainelll978(l (i.e. in 
photons cm~ 2 s _1 sr _1 ), 9 is the angle of the considered direction with respect to the normal to the surface. The UV 
opacity tiooo has been expressed as a function of the visual extinction A pp through the ratio of UV-to-visible opacities 
(tiooo/t550o)j where the subscript stands for the mean wavelength in A. 
The UV field at point r within a sphere with a r~ 2 density law is 



47r Jjjv (r) 



2ttGo I sin9exp{ 
'o 



T 100 o(r,9)}de 



T1000(V, 



'1000' 



il2 



atani 



vr 



(B2) 



where it is necessary a boundary conditio n, e.g. the UV optical depth t^qqq at point r* along the radial direction. 

We have assumed Tiooo( r =10 )=12.7 ijDotv fc Leund 11998) . although the opacity at UV wavel engths is highly 
uncertain for IRC+10216, and a (T1000/75500) ratio of 3.1 as measured by [Rouleau &: Martini l|1991f) for amorphous 
carbon grains of 0.05 /im, although different values up to 8 can be found in the literature depending on grain nature 
and size. The radial profile A v (r) ob tained is shown in F ig. IB1I where it is compared to the standard relation A v ccNh 
found for the interstellar medium bv lBohlin et alJ l)1978[L The interstellar UV field G has been decreased by a factor 2 
with respect to the standard Go to make the abundance distributions of C2H and CN (photodissociation products of 
C2H2 and HCN respectively) to peak at ~ 15" and ^20" respectively in agreement with intcrferometric observations 
ijDaval fc Bieginjll995l:lGuelin et al.1 1999). For our assumed distance of the star (150 pc) C2H peaks at 3.4xl0 16 cm 
and CN at 4.5x 10 16 cm (see Fig.|SJ|. We point out that this procedure is affected by a degeneracy between at least two 
parameters: the star distance, which transforms radial distances to angular distances, and the UV field and/or the A v 
radial profile, which make molecules photodissociate at a shorter or larger radial distance. Despite the uncertainties 
in the knowledge of these quantities, and even if our assumption of a less intense UV field than the standard is not 
correct, the abundances of the different molecules will be predicted at the right angular positions. 
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